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Abstract 



I In the study of many-particle systems both the interaction of particles can be 

essential and such feature as their internal (composite) structure. To describe these 
aspects, the theory of deformed oscillators is very efficient. Viewing the particles as 
p, g-deformed bosons, in the corresponding p, g-Bose gas model we obtain in explicit 
form virial expansion along with the 2nd to 5th virial coefficients. The obtained 
virial coefficients depend on the deformation parameters p, q in the form symmetric 



O I under p q, and at p — )■ 1, (7 — )■ 1 turn into those known for usual bosons. Besides 



real parameters, we analyze the case of complex mutually conjugate p and q and find 
interesting implications. Also, the critical temperature is derived (for the p, g-Bose 
gas) and compared with the Tc of standard case of bosons condensation. Similar 
results are presented for the deformed Bose gas model of the Tamm-Dancoff type. 



^ ' Keywords: deformed Bose gas model, virial expansion, virial coefficients, 

s 



critical temperature 

§■ 1 Introduction 



Diverse models of deformed Bose gas, elaborated on the base of a set of some deformed 



>• \ oscillators or deformed bosons, appeared in the literature from early nineties, see e.g. 

^ ■ Refs. |1]-|H]- Some analysis of a number of works on this topic has been done in Ref. |9]. 

^ ■ Since then, the statistical mechanics of a gas of q- or p, g-bosons was extensively explored 

^ . (and further extended) as witnessed by |10)-|19j and many others. 

With the aim of treating possible intermediate statistics behavior of physical system, 
one can use the method based on a deformation of quantum algebra of creation, anni- 
hilation and number operators. Deformed oscillator algebra is a generalization of the 
harmonic oscillator or Heisenberg algebra. Besides intermediate statistics, there exist 
^ . other convincing reasons of why one should use some model of deformed oscillator instead 

H ! of the usual quantum harmonic oscillator. The models of deformed oscillators are efficient 

in description of interacting particles system, or if the finite proper volume of particles, 
or their composite nature (internal structure) are to be taken into account. The theory 
of g-oscillators is related to the theory of quantum groups as originally shown in Ref. |20] 
and Ref. [21] . As a direct generalization of g-oscillators, in Ref. |22] and Ref. |23] the 
two-parameter family of p, g-oscillators has been introduced. 

The two-parameter p, g-deformed analog of Bose gas model can be developed, see 
[12] . |13| . using the set of p, g-deformed oscillators. Some further aspects of that model 
were studied in the works [11], [15]. In our present work we explore the thermodynamics 
of (some version of) p, g-Bose gas model. In order to obtain certain thermodynamical 
quantities of the deformed Bose gas we use the two-parametric generalization of g-calculus 
realized by replacing the ordinary derivative with the p, g-analog of Jackson derivative (the 
latter is known as g-deformed analog of usual derivative). 
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For the p, g-deformed Bose gas we investigate two regimes: high temperatures and 
small densities (A^/u ^ 1), or low temperatures and large densities (A^/w ^ 1). In 
the first case we obtain explicitly the virial expansion of the equation of state and find 
new virial coefficients in addition to few known ones (mainly in more special cases of 
one-parameter deformed Bose gas models [7], [IH])- The virial coefficients reflect effective 
inter-particle interactions and in the considered case depend explicitly on the deformation 
parameters p and q. From our formulas, the ordinary boson gas results can be recovered 
in the corresponding limit p — )■ 1, g — )■ 1. On the other hand, setting p and q as complex 
conjugates: p = q = re^^ , we gain yet another presentation of the virial coefficients. In that 
case we examine the behavior of the difference between deformed coefficients and their 
ordinary non-deformed prototypes (corresponding to non- deformed Bose gas) as function 
of the new parameter 6, at fixed values of the parameter r (the modulus). The situation for 
which Bose-Einstein condensation does occur in the considered two-parameter generalized 
boson gas is studied when X^/v ^ 1, and in this case the critical temperature Ti^'^^ of the 
condensation of deformed bosons is obtained. Comparing this with the of usual Bose 
gas we deduce that the critical temperature of the p, g-Bose gas is larger than the critical 
temperature of usual bosons, for the considered range of parameters. The dependence of 
the ratio Ti^''^^ Tc on the deformation parameters p, q is explicitly studied, visualized by 
the corresponding figures and the distinction from the results in Ref. [I5] indicated. We 
also consider similar results for the g-deformed Bose gas model of Tamm-Dancoff (TD) 
type, which constitutes a distinguished special case. 

2 Deformed oscillator models 

In the theory of deformed oscillators it is convenient to use the concept [TT] of deforma- 
tion structure function (p{N) which determines the particular deformed oscillator model 
described by its respective oscillator algebra. With (p{N) = a^a, the generating elements 
a, a^, of the algebra obey the relations 



[iV,at] 



[N,a] 



= —a, 



(1) 



aa'^-a^a = ip{N+l)-Lp{N), aa"^ = ip{N+l). 
In deformed Fock space the ground state vector obeys usual relations 



(2) 



a|0) = 0, iV|0)=0, (0|0) = 1 



(3) 



and for the n-particle excited state in this (y9(n)-extended Fock space we have 



N\n) = n\n), (p{N)\n) 



(4) 




(5) 



a^\n) = ^y ip{n+l)\n+l) , a\n) = ^/(p{n)\n — l) . 



(6) 
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In this paper we deal with the system of so-called p, g-bosons {p, g-deformed oscillators) 
for which the deformation structure function (f{n) in (12])- © is 

with [n\p^q being the p, g-number corresponding to a number n. Let us note that from the 
two-parameter family one infers a plenty of different one-parameter (g-deformed) models 
by imposing different relations p = f{q), see |26] for more details. 



3 Thermodynamics of the gas of deformed bosons 

In this Section we deal with a particular variant of the p, g-deformed Bose gas model. In 
the grand canonical ensemble, the Hamiltonian of the system is taken as 

H = Y,{e^-fi)N, (8) 

i 

where Si is the (kinetic) energy of particle in the state labeled by i, Ni is the boson number 
operator relative to Ei and /i is the chemical potential. One-particle non-relativistic ener- 
gies are given as Ei = p'j/2m where p = \']^\. Thermal averages can be calculated using 
the usual formulas of quantum statistical mechanics. Therefore, the thermal average of 
an operator A is given by the standard formula 

{A)=Tr{pA), P = ^ (9) 
where p is the equilibrium statistical operator, Z = Z{z) is the partition function, 

Z = Tr{e~^''). (10) 

Besides, (3 = 1/kT, T is the temperature; the Boltzmann constant k is put k = 1. 

To develop thermodynamics of the gas of deformed bosons we start with the logarithm 
of grand partition function 

InZ = -^ln(l-ze-'^"0 (11) 

i 

where z = e^^ is the fugacity. For obtaining thermodynamic functions of g-deformed Bose 
gas the usual derivative is replaced by the Jackson or g-derivative pi] 



^/(,)./MzM, (12) 

ax qx — x 

Note that in the limit g — 1 the usual derivative is recovered from the g-derivative. 

For the algebra given by ([T]), (l2]| and (f{n) from ([7]) we use the two-parameter p, g- 
generalization of Jackson derivative! being a direct p, g-extension of ( lT2l) . i.e. 



px — qx px — qx 



Let us note that in Refs. [H], [TS] yet another two-parameter analog D^''^^ of the Jackson derivative 

Inp— Inq ^ 



is exploited which is related with oi''''^'' in ^ as follows: D^^''^^ ^ . ^Zf J':^''^^ 
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When p = 1, the latter p, g-extension reduces to the Jackson derivative (fT2i) . 

Then, the total number of particles in the p, g-deformed theory is determined as 

jY = ^-p(p,g) i^z=^— ln{l-qze-^'^)-J2H'^-PZ(i~'^''} ■ (14) 

^ ^ i i 

For large volume V ^oo and large number of particles A^, we replace summation over 
level i by integration over the 3-momentum k. Note that the two sums in (HM diverge at 
z = 1 and p — > 1, g — )■ 1 when ei = 0. That is, these £i = terms can be as large as all the 
rest of terms in the sums. For that reason the = term is isolated, while the remaining 
two sums are replaced by integrals: 

Going over in k -space to spherical coordinates if,6,k = {k-k )^/^, the integration is 
performed with the result given through polylogarithms. Since the latter can be presented 
in the form of series, we derive the result 



-4 



p — q \-' — ' r 



2-^ ^5/2 2-^ ^5/2 



r=l 



p-q 



■pz) ^ 

r=l 



where no = ^;z^lii(Yr^) and A= y is the thermal wavelength. Now rewrite Eq. (1161) 
as 



1 1 ~ . ^ , "0 
- = 3^^73/2(^;p,g) + y 

where v= V / N , and the general formula for the (p, g)-extended function 'cjn reads: 

1 



(17) 



9niz;p,q)- 



p — q \ ^ — ' r 

^ r=l 



{zpY y. {Zqy \ _ I' Jp. 



r=l 



^ I T) 



+1 



(18) 



r=l 



This function is nothing but p, g-deformed generalization of the gn{^z) function known say 
from Ref. 1251: 



1=1 



(19) 



The function (ITS]) is real for real p and g. We assume 0<p<l,0<g<l. In what 
follows we will use 5(3/2 and 5^5/2. Namely, we will exploit (for small z) their power series 
expansions, say up to order. 

Then the thermodynamic relation PV/ T = InZ can be rewritten as 



j; = j^95/2{z;p,q) + -ln(l -2;). 



(20) 
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3.1 Virial expansion of the equation of state 

In the case of high temperatures and small densities {X^/v ^ 1), the average distance 
v^^^ between particles is much larger then the thermal wavelength A. In this case the 
quantum effects, the no, and the second term in (I20p can be viewed as negligibly small. 
Then, if the temperature of the gas of p, g-deformed bosons is high, so that T > tJ^'*^^ 
(see Sect. 1 below for ri^'*^^), from ([HD and ([20]) we find 

1 1 _ Pv V ^ 

- = p93/2{z;p,q), -Y=p95/2{z;p,q). (21) 



From the first equation in fl2T]) with account of ffTSl) we obtain 

^ 05/2 ^ Q5/2 ^ /I5/2 ^ c;5/2 ^ y"^"^' 



y 25/2 ' 35/2 ' 45/2 ' 55/2 

Inverting the latter equation yields 



V 



2^/^\v V 24 35/2 \ 215/2 25/235/2 25 



, mL , 7[2]L[3]p.g [3]L , 3[2],,J4],,, [5],,AA^V .,o^ 

29 2533/2 34 213/2 +53/2^1^^^- ^^-^^ 

From the second equation in ( l2T]l . with 515/2(2;; p, g) from f fTSl) taken in the form expanded 
in z, by the use of f l23l) we derive for the equation of state of the deformed p, g-Bose gas 
the desired virial expansion 

where the virial coefficients from 2nd to 5th read: 



A = -Mm 5 = Mm _ ^Bm c - 3[4],,, 5[2]3 



27/2 ' 25 37/2 ' 25/233/2 27 217/2 ' 

p mL [^]U^U . [2]p,j4],,, 2[3]^,^ 4[5],,, 
210 2^35/2 211/2 + 35 57/2 • 

Note, in the just obtained expressions for virial coefficients, p and q appear only through 
{p, g)-numbers [2]p^q, [3i]p,q, etc. Due to that, the exchange p q symmetry remains intact. 
Let us note that the above expressions for A, B, C, D yield the results for the particular 
cases of g-bosons. Namely, at p = 1 these reduce to the virial coefficients of the AC type 
g-Bose gas, and putting p = q~^ yields the BM type g-Bose gas virial expansions and virial 
coefficients. Another distinguished TD type case is considered separately in Sec. [5] below. 
Also, it should be noted that by putting p = f{q) where /(g) is some fixed function, as 
considered in Ref. [26], we readily obtain the relevant results (virial coefficients) for the 
corresponding, to this choice of /(g), version of g-deformed Bose gas model. 

Remark that the one-parameter (g-deformed) version of the above A, B, C was given 
in [18], [7]. The p, g-deformed 3-rd, 4-th and 5-th virial coefficients B, C and D are 
new. Note also that our A = A{p,q) differs from the respective coefficient in [14]. In the 
limiting no-deformation case p = q = 1 from the deformed virial coefficients A, B, C, D 
one recovers virial coefficients Aq, Bq, Cq, Dq of the usual Bose gas: 



25/2' 23/231/2 25 211/2' 
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35/2 ' u 26 22 ■ 33/2 25/2 33 55/2 • 

It is seen that Aq, Bq and Co coincide with the well-known virial coefficients of Bose gas 
given in textbooks [27], |28| . as it should. 



3.2 Complex parameters of deformation and virial coefficients 

Now take the deformation parameters q as complex, mutually conjugate ones: 

p = re'\ q = re"^^. (25) 

Then, instead of p and q we have the parameters r and 6 in terms of which the coefficients 
of virial expansion ( l24l) take the form: 

~ rcos9 ~ r^cos^B Ir^VlcoslQ + 1) 

A = TTT^, B ^ ^ 



25/2 ' 23 37/2 



~ r^cos9{2cos2e + 1) 'ir^^cosW + cosO) hr cos 6 

^ 23/233/2 ^6 2"/2 ' 

~ _ Ir^cos^e r^cos^9{2cos29 + 1) r^cos9{cos39 + cos9) 
- 22 . 35/2 + 2^2 

2r^(2cos2^ + 1)2 4r^(2cos4^ + 2cos2^ + 1) 
35 5^2 ■ 

Note, at r = 1 these formulas again reduce to respective results for the g-Bose gas of BM 
type g-bosons, with the phase like deformation parameter g = e*^. 

Physical meaning of the complex deformation parameters fl25|) can be commented as 
follows. If one deals with the usual non-ideal Bose gas, the virial coefficients contain terms 
responsible for the (two-particle, three-particle etc.) effective interactions. When dealing 
along the same lines with the gas of deformed bosons (deformed Bose gas), we gain that 
the substructure of particles or additional inter-particle interaction is effectively taken into 
account. Then the ffist non-trivial virial coefficient A or A reflects modified two-particle 
interaction, the second one B or B involves modified three-particle interaction and so 
on. With two parameters r and 9 at hands we might hope it is possible to diminish to 
zero two chosen types of interaction, say, the two-particle interaction together with the 
three-particle one (this would happen if A = and B = simultaneously). However, 
the parameter r (the modulus) turns out to be of no help for that aim. So, using the 
remaining parameter 9 we can find its value (s) for which only one of the two: A or B can 
be made zero. If that happens say for A, we may conclude: at this value of 9 we encounter 
mutual compensation of the two-particle interaction present in the usual non-ideal Bose 
gas against the additional contribution due to deformation, i.e. due to the physical reason 
(finite proper volume or substructure of particles) effectively taken into account by the 
deformation. 

On the other hand, we can find the measure of deviation of "deformed coefficients" 
from the known virial coefficients [27], |28) of the standard Bose gas: 

a = A-Ao, (3 = B-Bo, = C - Co, 6 = D - Dq. (26) 

In Fig. [Tjwe plot the dependence of a, /3, 7, 5 defined in (126|) on the parameter 9 at some 
fixed values of the second parameter r. 
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Figure 1: Dependence of a, /3, 7, 6 on the parameter 6* at r = 0.75, 1, 1.25. 

As seen, for different values of 6 the differences a, (3 and 7 can be positive, negative 
or zero, whereas the difference 5 is always negative (i.e., the net contribution to the 
effective five-particle interaction due to deformation is smaller than the interaction in the 
non-deformed Bose gas). 

In Fig. [2] we plot the dependence of a, /3, 7, 5 on the parameter 9 at fixed value r = 1 
(that implies the deformed model of BM type g-bosons). 



Diff.,A, 




Figure 2: Dependence of the differences Aj = (a, /3, 7, 5) from Eq. fl26|) on the parameter 
6* at r = 1. 

Remark. From FigslH [2] we observe remarkable thing: for all ^ 7^ 0, we have A > Aq, but 
D < Dq. That is the effect of deformation significantly increases attractive two-particle 
interaction for most of values of parameter ^ (at r < 1 and 6^0). This is in some 
accord with remark in [TU] for the case of quon-based deformation of Bose gas. In the 
domain of five-particle interaction the picture is opposite - the deformation modifies the 
interaction towards weakening. Although this regime (of high temperatures) is opposite 
to that one where the critical Tc is actual (see next Section), we nevertheless note that 
the deformation better promotes Bose condensation that results in increasing (due to 
the increase of attractive two-particle interaction). This fact may be of potential applied 
value. 
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4 Bose condensation in g-Bose gas and critical tem- 
peratures 

Here we study the opposite situation (to the Subsection 3.1), namely, the case of low 
temperatures and large densities [X^/v ^ 1). Eq. fl2T]) gives the equation of state for 
p, g-deformed Bose gas consisting of non-relativistic particles of mass m confined in the 
volume V. To explore the equation of state in detail we have to find the fugasity z as a 
function of temperature, and the specific volume should be found by solving the equation 
(fT7|) that involves besides no = ^^In^YZ^j) ^l^o v=V/N and the thermal wavelength 
\=^/2TxWJmT. Rewrite Eq. ^ as 

A^y = ^-^3/2(^;p,g). (27) 

We see that the value tiq/V is positive if the temperature and specific volume obey the 
following inequality: 

^> ^3/2(1; P,g). (28) 

In other words, the large though finite number of particles occupies the lowest energy level 
with £j = (ground state). That is, the phenomenon of Bose condensation takes place. 
For given specific volume v, the critical temperature can be found from the equality 

\l = vg^l2{l-,p,q). (29) 
Then, the critical temperature Tc^''^^ of the p, g-deformed Bose gas results as 

Moreover, we obtain the relation between the critical temperature of the considered p, q- 
deformed Bose gas and the Tc of usual gas of bosons in the form of the ratio: 



( 2.61 



> ^3/2(1; p,g), 



(31) 



Note that in the no- deformation limit p — )■ 1, g — )■ 1, the function (?3/2(l;P, q) in ( 130|) goes 
over into (73/2(1) = (^(|) = 2.61, as seen from ( !T9|l at 2; = 1. Due to this, the critical 
temperature J'i^''^'' of deformed Bose gas reduces to the critical temperature Tc of usual, 
non-deformed Bose gas and thus the ratio is (t!^''^^/ Tc)\p=q=i = 1. 

In Fig. ^{left) we give the plot of (13T]) as a function of the deformation parameters p, q 
such that p < 1, q < 1. This (convex upwards) behavior with respect to p and q differs 
from that in jT3]. Besides, this ratio increases with increasing (in the both parameters p, 
q) extent of deformation measured by 1 — p and 1 — q. In analogy with the above case 
of high temperature and low density, where we considered both real p, q and the variant 
with complex deformation parameters p = q = re^^, here in Fig. [3] (right) we also present 
the picture for (Tc^'^^ Tc) depending on complex p, q through the modulus r and the 
phase 9. 

As seen, the critical temperature Tc^'''^ for the p, g-deformed Bose gas is larger than 
the critical temperature Tc for the non-deformed boson gas, at least in the chosen region 
P < 1, Q' < 1 of the deformation parameters (see also Remark ending Sect. |3]). This fact 
may play an important role in future analysis involving real gases. Even more interesting 
is the picture for the ratio in Eq. fl3Tl) (at p = g = re*^ versus the parameters r and 6). 
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Figure 3: Left: The ratio Tc /Tc of the critical temperatures given by Eq. (1311) as a 
function of the deformation parameters p, g such that < p < 1 and < g < 1. Right 
The ratio T^^'^^ versus the deformation parameters r and 6', 0<r<l, 0<6'< 2tx. 



5 Deformed Bose gas of Tamm-Dancoff type 

Let us consider thermodynamics of the g-Bose gas of Tamm-Dancoff (TD) type. Recall 
that the defining (along with ([1])) commutation relation 

aa) — qa)a = (32) 

for the TD g-bosons, and the g-bracket of TD type 

[X]td = Xg^-i (33) 

stem from those of the two-parameter g-deformed model if one puts p = q, see ([7]). We 
assume that < g < 1 in the case of TD bosons. For some exotic properties of the TD 
type deformed oscillator see Ref. |29] . 

Performing when X^/v ^ 1 same calculations as in Sec. [3l in the Tamm-Dancoff case 
we find the coefficients for the virial expansion 



^™ - - ^^0' ^™ - ^ (^23/231/2 " 25 ~ 21172 J - ^ '^O' 

BtD = ^'(^ - ^) = ^'^0, ^TD = - 22 .33/2 + ^ + § - ^) = ^'^0- 

It is worth to note that, due to g < 1, the TD Bose gas virial coefficients are lowered with 
respect to non-deformed virial coefficients, by the factor of g'^~i in the /c-th virial coef- 
ficient. Accordingly, this weakens in (^~^ times the usual /c-particle effective interaction 
within the TD g-Bose gas. 

In the case of low temperatures and large densities, }? jv ^ 1, the critical temperature 
of Bose gas of Tamm-Dancoff type is 



T™ / 2.61 



■ c 



(34) 



where g^j^^{\\ q) is the TD analog of function ( JTOl) at z = l and n = 3/2. In general, 

00 r— 1 r 

= (36) 

r=l 
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Figure 4: The ratio Tj^ jT^ of the TZ^-type deformed critical temperature Tj^ to the 
non-deformed one Tc versus the deformation parameter q where < g < 1 . 

In Fig. m we plot the ratio versus deformation parameter g for < g < 1. As seen, 
the larger is deformation (deviation 1 — g from the non-deformed value g = l), the larger 
is the critical temperature T]:^ of deformed bosons of TD type. At g = 1 the ratio (15^ 
gives 1, as it should be. 

6 Concluding remarks 

In this paper we developed some version of the p, g-deformed analog of the Bose gas 
model (p, g-Bose gas model). We studied the thermodynamics of such a gas at high (low) 
temperatures and low (large) densities. To obtain the thermodynamical quantities we 
utilized the extension of g-calculus based on the direct p, g-generalization of the Jackson 
derivative. 

For high temperatures [}? j v <^ 1), dealing with the equation of state we have obtained 
the virial expansion and correspondingly derived the nontrivial virial coefficients A, 5, C, 
D all depending on the deformation parameters p and g. Note that these parameters are 
contained in the virial coefficients through the p, g-numbers only and are thus symmetric 
under the exchange p -f-)- g. In the limit p = g = 1 the virial coefficients Aq, -Bqj C*o, -Dq of 
usual Bose gas are recovered. The parameters p and g are taken first to be real, and then 
also as complex valued such that p = g = re*^. In these new variables r and Q we analyzed 
the virial coefficients as well. The differences between deformed coefficients and their 
non-deformed counterparts have been studied from the viewpoint of their dependence 
on the parameter Q. Figures [H |2] demonstrate: there exist some special value(s) of the 
6'-parameter (and r) for which any chosen difference from f l26l) can vanish due to some 
mutual compensation. On the other hand, for special fixed value Q (at certain r), the 
virial coefficient A, and thus the effective two-particle interaction can be vanishing. Note 
that the effective two-particle interaction (and vanishing of it) means superimposing (and 
mutual compensation) of the conventional two-particle interaction, present for the usual 
non-deformed non-ideal Bose-gas, and the additional interaction imported due to the 
deformation used by us. Evidently, the same can be said about each from the rest of 
virial coefficients (taken alone). 

At low temperatures (A^/ w ^ 1) we find both the critical temperature t'^'^^ of deformed 
p, g-Bose gas and the explicit dependence of the ratio T^f^'^^j Tc on the parameters p and 
g (either taken to be real, or as the pair of mutually conjugate complex values). Similar 
results are deduced for the deformed Bose gas of Tamm-Dancoff type. We have found, for 
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the whole range of the considered p, q, that Ti^'*^^ >Tc (such inequahty retains in the case 
of TD-type g-Bose gas also). We hope this feature may have interesting consequences. To 
conclude, the issue of Bose condensation in its p, g-deformed manifestation (with Ti^'*^^ or 
Tc^'^^ higher than T^) is reproduced and pictured in Fig. 3, left and right. It would be 
interesting to make comparison of the obtained results with the critical Tc of a real gas. 

Acknowledgments 

This research was partially supported by the Special Program of the Division of Physics 
and Astronomy of the NAS of Ukraine. 

References 

[I] M. Martm-Delgado, J. Phys. A24, (1991) L1285. 

[2] R. Bullough, N. Bogoliubov, Phys. Lett A168, (1991) 801. 

[3] G. Agarwal, S. Chaturvedi, Mod. Phys. Lett. A7, (1992) 2407. 

[4] T. Altherr, T. Grandou, Nucl. Phys. B402, (1993) 195. 

[5] P. Neskovic, B. Urosevic, Int. J. Mod. Phys. A7, (1992) 3379. 

[6] M. Chaichian, R. Gonzalez Felipe, C. Montonen, J. Phys. A26, (1993) 4017. 

[7] M.R.-Monteiro, I. Roditi, L.C.M.S. Rodrigues, Mod. Phys. Lett B7, (1993) 1897. 

[8] Gang Su, Mo-lin Ge, Phys. Lett A173, (1993) 17. 

[9] S. Vokos, C. Zachos, Mod. Phys. Lett. A9, (1994) 1-10. 

[10] S.S. Avancini, G. Krein, J. Phys. A: Moth. Gen. 28, (1995) 685. 

[II] S. Meljanac, M. Milekovic and S. Pallua,Phys. Lett. B 328, 55 (1994). V.I. Man'ko, 
G. Marmo and F. Zaccaria, Phys. Lett. A197, (1995) 95. D. Bonatsos and C. Daskaloy- 
annis. Prog. Part. Nucl. Phys. 43, 537 (1999). 

[12] L.V. Adamska, A.M. Gavrilik, J. Phys. A37, (2004) 4787. 

[13] A.M. Gavrilik, SIGMA 2, (2006) 074. 

[14] A. Algin, Commun Nonlinear Sci Numer Simulat 15, (2010) 1372. 

[15] A. Algin, J. Stat. Mech., (2008) P10009. 

[16] A. Lavagno, P. Narayana Swamy, Phys. Rev. E65, (2002) 036101. 

[17] M. Rego-Monteiro, L.M.C.S. Rodrigues, S. Wulck, Phystca A259, (1998) 245. 

[18] A.M. Scarfone, P. Narayana Swamy, J. Phys. A41, (2008) 275211. 

[19] A.M. Gavrilik, A.P. Rebesh, Eur. Phys. J. A47, (2011) 55. 

[20] L. Biedenharn, J. Phys. A22, (1989) L873. 



12 

[21] A. Macfarlane, J. Phys. A22, (1989) 4581. 

[22] A. Chakrabarti and R. Jagannathan, J. Phys. A: Math. Gen. 24, (1991) 711. 
[23] M. Arik et ai, Z. Phys. C55, (1992) 89. 
[24] F. Jackson, Mess. Math. 38, 57 (1909). 

[25] K. Huang, Statistical mechanics (John Wiley & Sons, Inc., New York - London, 
1963). 

[26] A.M. Gavrilik, A.P. Rebesh, Mod. Phys. Lett. A, Volume 23, Issue 13, (2008) 921. 

[27] R. Kubo, Statistical mechanics (North-Holland Publishing Company-Amsterdam, 
1965). 

[28] R.K. Pathria, Statistical mechanics, (Second edition, Butterworth Heinemann, 1996). 
[29] A.M. Gavrilik, A.P. Rebesh, Mod. Phys. Lett. A 27, (2007) 949. 



